THE MAIN OBJECTIVE OF THIS PAPER IS TO PROVIDE AN EFFICIENT NUMERICAL METHOD FOR SOLVING THE ADVECTION EQUATION USED IN MIXTURE MULTIPHASE MODEL, TWO-DIMENSIONAL TIME SPLITTING METHOD. FOR ADVECTION DOMINATED PROBLEMS, SUCH AS THE MULTIPHASE FLOW THE EFFECT OF THE NUMERICAL DIFFUSION IN THE ADVECTION SCHEMES WILL ADVERSELY AFFECT THE ACCURACY OF THE NUMERICAL RESULTS. TO THIS END, FIRST SCHEMES WITH DIFFERENT ORDERS WERE EVALUATED. USING HIGHER ORDER SCHEMES INTRODUCED A REDUCTION IN NUMERICAL DIFFUSION. HOWEVER, SUCH SCHEMES LED TO UNREAL, OUT OF RANGE AND NON-MONOTONIC RESULTS. SUBSEQUENTLY, THE SOLUTIONS OFFERED FOR SOLVING THESE PROBLEMS BY OTHER RESEARCHERS IN THE PREVIOUS STUDIES WERE EXAMINED. WITH DUE REGARD OF THE MIXTURE MODEL CHARACTERISTICS, THE MOST EFFECTIVE OF THESE SOLUTIONS WERE THEN SELECTED AND SUBSEQUENTLY PRESENTED FOR SOLVING THE ADVECTION EQUATION IN THE MIXTURE MODEL.